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Abstract. Combining pointwise Green's function bounds obtained in a companion 
paper [MZ.2] with earlier, spectral stability results obtained in [HuZ], we establish 
nonlinear orbital stability of small amplitude viscous shock profiles for the class of 
dissipative symmetric hyperbolic-parabolic systems identified by Kawashima [Kaw], 
notably including compressible Navier-Stokes equations and the equations of magne- 
tohydrodynamics, obtaining sharp rates of decay in L p with respect to small L 1 f)H 3 
perturbations, 2 < p < oo. Our analysis follows the approach introduced in [MZ.l] 
to treat stability of relaxation profiles. 



Section 1. Introduction 

Consider the class of degenerate parabolic conservation laws of dissipative, sym- 
metric hyperbolic-parabolic type in the sense of Kawashima [Kaw], i.e., systems 

(1.1) G(U) t + F(U) x = (B(U)U x ) x , 
satisfying 

(Al) dF, dG, B symmetric, dG > 0, B > (symmetric hyperbolic-parabolicity), 

(A2) No eigenvector of dFdG~ x lies in the kernel of BdG~ x (dissipativity), 
and 

(A3) U=(™\ B= (® °\ b>0 (block structure), 

in some neighborhood U of a particular base point [/*, where U G IR n , u G M n_r , 
v G W, and b G IR rXr . As discussed in [Kaw], this class of equations includes 
many physical models arising in continuum mechanics: in particular, compressible 
Navier-Stokes equations and the equations of compressible magnetohydrodynamics 
(MHD). We make also the regularity assumption 

(HO) F, G, B G C 3 , 

needed for our later analysis. 

An interesting class of solutions of (1.1) are viscous shock profiles, or asymptoti- 
cally constant traveling wave solutions 

(1.2) U = U(x-st); lim U(z) = U± 

z— >±oo 
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connecting endstates U± corresponding to (discontinuous) shock waves of the asso- 
ciated hyperbolic system 

(1.3) G(U) t + F(U) x = 0. 

Existence of viscous profiles has been treated in the large for compressible Navier- 
Stokes by, e.g., Gilbarg [Gi] and for MHD by Conley and Smoller [CS] (see [Fre.l] 
for a treatment in the small amplitude case). Most recently, Freistiihler [Fre.2] has 
carried out a general treatment of existence in the small, analogous to that of Majda 
and Pego in the strictly parabolic case [MP], for arbitrary models (1.1) satisfying 
(A1)-(A3). 
Let 

a\(u) < ■ ■ ■ < a n (u) 

denote the eigenvalues of Aq := dFdG -1 ^), rj(U) and lj(U) a smooth choice of 
associated right and left eigenvectors, lj -rk = S 3 k , and assume within neighborhood 
U of base point U* that: 

(HI) The pth characteristic field is of multiplicity one, i.e. a p {U) is a simple 
eigenvalue of Ag{U). 

(H2) The pth characteristic field is genuinely nonlinear, i.e. Va p • r p (U) ^ 0. 

(H3) The eigenvalues of dFudG^ (necessarily real and semisimple) are: (i) of 
constant multiplicity; and, (ii) different from a p . 

Then, we have: 

Proposition 1.1 ([Fre.2]). Let (Al)-(AS), and (Hl)-(HS) hold, and F, G, B e 
C 2 in (1.1) (implied in particular by (HO)). Then, for left and right states U± lying 
within a sufficently small neighborhood V C U of U*, and speeds s lying within a 
sufficiently small neighborhood of a p (U*), there exists a viscous profile (1.2) that is 
"local" in the sense that the image of u(-) lies entirely within V if and only if the 
triple (U-,U + ,s) satisfies both the Rankine-Hugoniot relations: 

(RH) s[G] = [F], 

and the Lax characteristic conditions for a p-shock: 

(L) a p (U-) > s > a p (U+); sgn (a,j(U-) - s) = sgn (aj(U+) - s) ^ for j ^ p. 

(Note: The structure theorem of Lax [La,Sm] implies that (RH), always a necessary 
condition for existence of profiles, holds for U± G V only if s lies near some a,j(U*); 
thus, the restriction on speed s is only the assumption that the triple (U-, U+,s) be 
associated with the pth and not some other characteristic field). 

Remark 1.2. The generically satisfied conditions (HI) and (H2) are implied by 
strict hyperbolicity and genuine nonlinearity, respectively, of the associated hyper- 
bolic system (1.3). In particular, they hold always for the equations of compressible 



CORRADO MASCIA AND KEVIN ZUMBRUN 



3 



gas dynamics, under the assumption of an ideal gas [Sm]. Condition (H3)(ii) corre- 
sponds to the requirement that "hyperbolic," or unsmoothed modes in the solution 
be noncharacteristic, at least as regards the principal characteristic speed a p . This 
technical condition is needed in order that the traveling wave ODE be of nonde- 
generate type, and is a standard assumption in the theory. See assumption (4) 
of [Fre.2], assumption (HI) (ii) of [Z.3], Appendix A. 2, or assumption (H5) of [SZ] 
for restatements in various different contexts. Condition (H3)(i) is an additional 
technical assumption that was used in the detailed Green's function analysis carried 
out in [MZ.2] for slightly more general systems; we suspect that it can be dropped 
in the present, symmetrizable case. Conditions (H3)(i)-(ii) are likewise satisfied for 
compressible ideal gas dynamics, for which dFndG±± is 1 x 1 and identically equal 
to particle velocity (u in Eulerian coordinates, in Lagrangian coordinates). 

Stability of viscous profiles has been examined for compressible Navier-Stokes 
equations in [MN,KMN,L.2], vith various partial results concerning special (mainly 
zero-mass) initial data; the results of [KMN,L.2] are restricted to small amplitude 
profiles, while the results of [MN] in the case of an isentropic 7-law gas apply to 
profiles of amplitude 01(7), with a — > 00 as 7 — > 1. Yet, as pointed out recently in the 
fluid-dynamical survey [Te] , the basic problem of stability of gas-dynamical shocks 
in its full generality remains open even for small amplitude waves, a significant 
gap in the theory of compressible gas dynamics. More recently, Humpherys and 
Zumbrun [HuZ] have established strong spectral stability of general, small amplitude 
shock waves of Kawashima class systems, of the type constructed by Freistiihler, 
generalizing a corresponding result of Goodman [Go. 1-2] in the strictly parabolic 
case (see Section 2, below). As discussed in [ZH,HuZ], strong spectral stability is 
roughly equivalent to, but slightly weaker than stability with respect to zero-mass 
perturbations. 

In [MZ.2], applying the general machinery developed in [ZH,MZ.l], we have 
shown, for a slightly more general class of systems and for profiles of arbitrary 
amplitude and type, that strong spectral stability, plus hyperbolic stability of the 
corresponding ideal shock (always satisfied for weak Lax shocks satisfying (Hl)- 
(H3)), are necessary and sufficient conditions for linearized orbital stability, and, 
moreover, yield extremely detailed pointwise bounds on the Green's function of the 
linearized evolution equations, analogous to those obtained for relaxation profiles in 
[MZ.l]. In combination with the result of [HuZ], these results imply L 1 n L p — > L p 
linearized orbital stability of small-amplitude shock profiles of Kawashima class sys- 
tems, with sharp rates of decay for all 1 < p < 00. 

The purpose of the present paper, extending the work of [HuZ, MZ.2], is to es- 
tablish L 1 n H 3 — > L p nonlinear orbital stability as solutions of (1.1) of small am- 
plitude profiles of general Kawashima class systems, with sharp rates of decay for 
all 1 < p < 00. More precisely, we shall establish: 

Theorem 1.3. Let there hold (Al)-(A3) and (H0)-(H3) for a general relaxation 
model (1.1), withU* andU as above. Then, forV C U sufficiently small, the viscous 
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profiles U described in Proposition 1.1 are nonlinearly orbitally stable from L 1 OH 3 
to L p , for all p > 2. More precisely, for initial perturbations U$ := Uq — U such 
that \Uo\Li n H3 < (o, Co sufficiently small, the solution U = (u,v)(x,t) of (1.1) with 
initial data Uq satisfies 

(1.4) \U(x,t) - U(x - S(t))\ LP < CCo(l + f)-** 1 " 1 /*) 
and 

(1.5) | (v(x,t) - v(x - S(t))) x \ LP < CCo(l + t)-^ (1 " 1/p) , 
for all 2 < p < oo, for some 5(t) satisfying 

(1.6) i<j(*)i<cco(i+*r* 

and 

(1-7) \S(t)\ <CCo- 

In particular, this implies nonlinear stability of all small amplitude gas-dynamical 
profiles under the assumptions of an ideal gas; see Remark 1.2 above. 

Our analysis in this paper, and in the companion paper [MZ.2], follows closely 
to that used in [MZ.l] to treat stability of relaxation profiles, making use of struc- 
tural similarities between the two problems pointed out in [Ze.2,Z.3]; in turn, the 
argument of [MZ.l] makes substantial use of the pointwise semigroup machinery 
introduced in [ZH,Z.2] to treat the parabolic case. As are general features of the 
pointwise semigroup approach, we obtain through this program extremely detailed 
pointwise bounds on the Green's function (more properly, distribution) of the lin- 
earized operator about the wave, sufficiently strong that the nonlinear stability 
analysis becomes comparatively simple. Moreover, these bounds (though not the 
nonlinear stability argument in which they are used, see Remark 1.5 just below) 
depend only on a generalized spectral stability (i.e., Evans function) condition and 
not on the amplitude or type (Lax, under-, or overcompressive) of the wave. 

Remark 1.4. The assumption of genuine nonlinearity (H2) is not needed either 
for the existence or the stability result, but is made only to simplify the discussion. 
Though we stated above only the restriction to the genuinely nonlinear case, exis- 
tence was in fact treated for the general (nongenuinely nonlinear) case in [Fre.2], 
substituting for the Lax entropy condition (L) the strict Liu entropy inequality (E) 
of [L.3]. Likewise, as described in [HuZ], the result of spectral stability may be ex- 
tended to the general case by substituting for the "Goodman-type" weighted energy 
estimates of [HuZ], Section 5, the variation introduced by Fries [Fri.1-2] to treat 
the nongenuinely nonlinear case for strictly parabolic viscosities. Accordingly, the 
result of Theorem 1.3 holds also in this case. 
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Remark 1.5. The restriction to small-amplitude shocks arises only through 
the energy estimates used to close the nonlinear iteration argument. In particular, 
it should be possible for the isentropic, 7-law gas case to treat amplitudes of the 
same order ot(j) treated in [MN]. A fundamental open problem is to remove the 
amplitude restriction altogether, as done in the strictly parabolic case [ZH,Z.2] 
and in the case of discrete kinetic relaxation models [MZ.l], replacing it with the 
generalized spectral condition of the linearized theory. 

Remark 1.6. Useful necessary conditions for viscous stability have been ob- 
tained in [Z.3] for arbitrary amplitude profiles of the more general class of models 
considered in [MZ.2], using the stability index of [GZ,BSZ]. Strengthened versions 
of the one-dimensional inviscid stability criteria of Erpenbeck-Majda [Er,M.l-3], 
these readily yield examples of unstable large-amplitude profiles, similarly as in the 
strictly parabolic case (see, e.g., [GZ,FreZ,ZS,Z.3]). This shows that the spectral 
stability requirement is not vacuous in the large amplitude case. 

Plan of the paper. In Section 2 we cite the spectral stability result of [HuZ], 
and in Section 3 the pointwise Green's function bounds of [MZ.2]. As an imme- 
diate corollary, we establish in Section 4 the linearized orbital stability of general 
shock profiles satisfying the necessary conditions of spectral and hyperbolic stabil- 
ity: in particular, of small amplitude profiles of Kawashima class systems. Finally, 
we establish in Section 5 the main result of nonlinear stability of weak profiles of 
Kawashima class systems, by a modified version of the argument of [MZ.l]. 

Note: Liu and Zeng have informed us [LZe.2] that they also have obtained 
nonlinear stability of weak Navier-Stokes profiles, by a different argument based 
on the approximate Green's function approach of [L.2]. This method is inherently 
limited to weak shocks of classical, Lax type, both at the linear and nonlinear level. 

Section 2. Spectral stability. 

Take without loss of generality s = 0, so that U = U(x) becomes a stationary 
solution. Then, the linearized equations of (1.1) about U take the form 

(2.1) U t = LU := -{A Q )-\AU) X + (A°)-\BU X ) X1 
where 

(2.2) B := B(U), A := dG(U), Av := dF(U)v - (dBv)U x ., 

Definition 2.1. We call the profile U(-) strongly spectrally stable if the linearized 
operator L about the wave has no spectrum in the closed unstable complex half- 
plane {A : Re A > 0} except at the origin, A = 0. (Recall, [Sat], that A = is 
always in the spectrum of L, since LU X = by direct calculation/differentiation of 
the traveling wave ODE) . 
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The spectral stability of small-amplitude viscous profiles of the slightly more 
more general class of dissipative, symmetrizable degenerate parabolic systems has 
been investigated in [HuZ]. To apply these results, we have only to note that (1.1), 
expressed with respect to variable G := G(U), becomes 

(2.3) G t + F(U(G)) X = (B(U(G))U(G) X ) X , 
or, in quasilinear form, 

(2.4) G t + F G^G X = {BG^G X ) X . 

System (2.4) is evidently symmetrizable (by the symmetric positive definite G^ 1 ), 
and dissipative (since the dissipativity condition is independent of coordinate sys- 
tem), with the additional block structure property that the left kernel of the new 
viscosity matrix BG^. 1 is constant. This is almost the class of equations considered 
in [HuZ], the difference being that there the block structure assumption was that 
the viscosity matrix have constant right instead of left kernel. This discrepancy is 
unimportant in the analysis, since, at the linearized level, one case may be con- 
verted to the other by the (linear) change of coordinates A°V := U, corresponding 
to similarity transform M — > (A°)~ 1 MA°. Alternatively, in the particular case of 
our interest, we may simply carry out energy estimates in the natural coordinates 
of (2.1); for related calculations, see Section 4.1, below. 

Thus, the results of [HuZ] apply to the somewhat larger class of dissipative 
symmetrizable systems of form (2.3)-(2.3), with the block structure condition that 
either the left or right kernel of the viscosity matrix be constant, and we may 
conclude, in particular: 

Theorem 2.2[HuZ]. Let (Al)-(A3) and (HO)-(HS) hold, and let U(x - st) be 

a viscous shock solution of (1.1) such that the profile {u(z)} lies entirely within a 
sufficently small neighborhood V CW ofU*, and the speed s lies within a sufficiently 
small neighborhood of a p (U*). Then, U is strongly spectrally stable, in the sense of 
Definition 2.1 above. 

Section 3. Pointwise Green's function bounds. 

In a companion paper to this one [MZ.l], we have investigated linearized stability 
and behavior of viscous profiles of systems of the general form (2.4), not necessarily 
symmetrizable, satisfying the standard set of conditions identified in [Z.3]. As 
remarked in [Z.3, MZ.l], these hold always for symmetric systems (1.1) satisfying 
(Al)-(A3), (H0)-(H3), and the small-amplitude profiles described in Proposition 
1.1; however, they may also hold in much greater generality, in particular for shocks 
of large amplitude, or nonclassical type. For simplicity of exposition, we shall 
restrict our discussion here to the present case of interest; for the general case, see 
[MZ.l]. 
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The first main result of [MZ.2], generalizing the corresponding results established 
for viscous, strictly parabolic, shocks in [ZH], and for relaxation shocks in [MZ.l], 
is: 

Theorem 3.1 [MZ.2]. Under assumptions (A1)-(A3), (H0)-(H3), small ampli- 
tude shock profiles are L 1 fl LP — > LP linearly orbitally stable for p > 1 if and only 
if they are strongly spectrally stable, with sharp decay bounds 

(3.1) \U(-,t) + 8(t)U'(-)\ LP KCil + ty^-^QUoUi + \U \ LP ) 

for initial data Uq G L 1 fl LP , some choice of 6. 

Theorem 3.1 is obtained (see, e.g., Corollary 4.2, below for one direction) as a 
consequence of detailed, pointwise bounds on the Green's function (more properly 
speaking, distribution) Q(x,t;y) of the linearized evolution equations 

(3.2) G t = LqG := —{AqG) x + (B G G X ) X , 

(3.3) Bq := B(A°)~ 1 (U(x)), Aqv := dF(A°)~ 1 v — dB(U(x))vU'(x), 

corresponding to (2.4), i.e., expressed with respect to the conservative variable G. 
We now describe these bounds, for use in the following sections. 

Let a^, j = 1, . . . (n) denote the eigenvalues of Aq(±oo), and Pf and as- 
sociated left and right eigenvectors, respectively, normalized so that lf l r^ = 5 J k . 
In case Aq is strictly hyperbolic, these are uniquely defined. In the general case, 
we require further that be left and right eigenvectors also of Pj°BqP^, 

Pp := R^L^ 1 , where L^ and denote mf x left and right eigenb locks as- 
sociated with the m^-fold eigenvalue af, normalized so that L^R^ = I m ±. (Note: 

The matrix PpB^Pp ~ L^BqR^ is necessarily diagonalizable, by simultaneous 
symmetrizability of Aq, Bq)- 

Eigenvalues a,j(x), and eigenvectors lj, rj correspond to large-time convection 
rates and modes of propagation of the degenerate model (2.4). Likewise, let a*(x), 
j = 1, . . . , (n — r) denote the eigenvalues of 

Aq := Aq, u - A G , 12 Bq^Bg, 21 
= A 11 {A° 1 )-\ 

and lj(x), r*(x) G R n_r associated left and right eigenvectors, normalized so that 

lj*rj = S J k . More generally, for an m*-fold eigenvalue, we choose (n — r) x m* blocks 
L* and R* of eigenvectors satisfying the dynamical normalization 

Lfd x R* = 0, 



(3.4) 
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along with the usual static normalization L* l Rj = S J k I m * ; as shown in Lemma 4.9, 
[MZ.l], this may always be achieved with bounded L*, R*. Associated with L*, R*, 
define extended, n x m* blocks 

(3 ' 5) ('•)• *!"(-b?X»*; 

Eigenvalues a* and eigenmodes £*, 1Z* correspond, respectively, to short-time hy- 
perbolic characteristic speeds and modes of propagation for the reduced, hyperbolic 
part of degenerate system (2.4). Note that our discussion in the introduction of 
condition (H3) is validated by the second equality in (3.4), a nontrivial consequence 
of symmetry/block structure in the original system (1.1). 
Define time-asymptotic, scalar diffusion rates 

(3.6) pf:=(l t j B G r j ) ± , j = l,...,n, 
and local, rrij x mj dissipation coefficients 

(3.7) r ] *{x):=-lfD*r*{x) : j = 1, . . . , J < n - r, 
where 

(3.8) D* := Ac, 12 Aq,21 - Ag,22Bq^ 22 Bg,21 + A g Bq^ 22 Bg, 21 

is an effective dissipation precisely analogous to the effective diffusion predicted by 
formal, Chapman-Enskog expansion in the (dual) relaxation case. As described 
in Appendix A2 of [MZ.2], these quantities arise in a natural way, through Taylor 
expansion of the (frozen-coefficient) Fourier symbol 

(3.9) -itA G (x)-eB G (x) 

of the linearized operator L G about £ = and £ = oo, respectively. 

The important quantities j3j were identified by Zeng [Ze.l,LZe.l] in her study 
by Fourier transform techniques of decay to constant (necessarily equilibrium) so- 
lutions (u, v) = (u±, v±) of relaxation systems, corresponding at the linearized level 
to the study of the limiting equations 

(3.10) U t = L G U := -A G U X + B%U XX 

as x — > ±oo of the linearized evolution equations (3.2). As a consequence of dissi- 
pativity, (A2), we have (see, e.g., [Kaw,LZe.l]) that 

(3.11) /3f>0, j = l,...,n 
and 



(3.12) Re cr(r]*{x)) > 0, j = 1, . . . , J < n - r. 
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Proposition 3.2[MZ.2]. For weak shock profiles, under assumptions (Al)-(A3), 
(H0)-(H3), the Green's function Q(x,t;y) associated with the linearized evolution 
equations (3.2) may in the Lax or over compressive case be decomposed as 

(3.13) Q(x, t; y) = H + E + S + R, 

where, for y < 0: 

J 

H(x,t;y) := Y^K*{x)q{y,t)5 x _- a * t {-y)Cf{y) 

(3.14) 

= j^K*{x)0{e-^)8 x -- ap {-y)Cf{y), 

3 = 1 

(3.15) 

E{x,t;y):= £ [4 t _]G'{x)l? ( errfn [ -±M) - errfn (-^M 

and 
(3.16) 



S(x,t;y) :=X{t>i} E r k l k {^^)~ 1/2 e 

a"<0 

+ X«>1> E r t -I t -'(*r^t)-'/> e -(-»-;«)VW« {-^\ 

a k >0 

+ *(<>.) E [^7t^^«)-V» e -<-«i>'/«i« (-j^lj) 

a fe >0, a j <0 



a fe >0,a+>0 



denote hyperbolic, excited, and scattering terms, respectively, and R denotes a faster 

decaying residual, satisfying: 

(3.17) 

R(x,t; y) = 0{e-^ x - y \ +t) ) 

n 

+ J20((t + l)"l/2 e -^ + + e -"l-l) t -V2 e -(x-y-a-t) 2 /Mt 
fc=l 

+ E X { |a- t |>H } ^((^l)- 1/2 ^ 1/2 )e- ( ^" (tH ^" l))2/Mt e^ 

a~>0, aj<0 



a fe >0,a+>0 
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(3.18) 

J 

R y (x, t;y) = J2 0(e-"%- a ;t(-y) + 0{e-^*-* + *>) 

3 = 1 
n 

+ O ((* + l)" 1 /2 e -^ + + e ~v\x\^ t -l e -{ X -y-a~ k tflMt 
fc=l 

a~>0,o~<0 

+ E X { |a fc -*|>| v | } 0((t+l)- 1/ V 1 )e-^-<C*Hv/- fc -|))V^ e -,*- j 

a~>0,at>0 

(3.19) 

J 

R x (x, t;y) = J2 0(e-^)6 x - ap (-y) + 0( e ->(l*-i'l+*)) 

3 = 1 
n 

+ J2<D((t + l)- 1 e-" a,+ + e-^l) t" 1 ^ + i)i/2 e -(x-y-« fc -t) 2 /Mt 
fc=i 

a">0, a7<0 

+ E X { | a - t |>| tf | } 0(t+l)- 1 /V 1 )e-(*--i + ('H^I)) a /^e-^-, 

a~>0,a+>0 

and 
(3.20) 



i2x»(x,t;y) = (d/3y)(^0( e ^%_ s * t (-y)) 

i=i 

J 

+ ^ 0(e-"%_ ffl ;t(-y) + 0(e-"f a! -w' + *>) 

3 = 1 
n 

+ J20((t + l)"3/2 e -^ + + e -*?M) t "3/2 (t + 1)e -( K - y - a -t) 2 /Mt 



fc=l 



+ £ x { |a fc -t|>|y|}0(t+l)- 1/2 t- 3 / 2 )e-^- o i"< t -l»/°*"l)) a / Mt e-'' a!+ 

a~>0,a7<0 

+ E x Ua ~ tl > lyl} o(t+ir 1/2 t- 3/2 )e-^^ 

a~>0,a+>0 

/or some rj, M > 0, where x ± denotes the positive/negative part of x, and indicator 
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function X{| a -t|>|j,|} *s one f or \ a kA — \v\ an< ^ zero otherwise. Symmetric bounds 
hold for y > 0. 

Here, the averaged convection rates a*Ay,t) in (3.14) an d (3.19)-(3.20) denote 
the time-averages over [0,t] ofa*(x) along characteristic paths z* = z*(y,t) defined 
by 

(3.21) dz*/dt = a*(z*), z*(0) = y, 

and the dissipation matrix Q = (j(y, t) G R m J' xm -? is defined by the dissipative flow 

(3-22) c*C*M = -r/*(**)C*, (*(y) = Im r 

Similarly, in (3.16), 

(3-23) z%^ := af (t - 



\ a k 



and 

(3.24) pf k (x,t;y):= 




represent, respectively, approximate scattered characteristic paths and the time- 
averaged diffusion rates along those paths. In all equations, aj, a*^ , lj, C*^ , 

Tj, TZf, /3 f and n* are as defined just above, and scattering coefficients [c% l _], 
% = — , 0, +, are constants, uniquely determined by 

(3.25) K-]r- + + [<-] " = r k 

aj<0 a+>0 

for each k = 1, . . .n, and satisfying 

(3-26) £ [cl_]r k = 14M = *, 

a~>0 a+<0 

where the constant vector tv is the left zero effective eigenfunction of Lq associated 
with the right eigenfunction G' . 

Proposition 3.2, the variable-coefficient generalization of the constant-coefficient 
results of [Ze.l,LZe.l], was established in [MZ.2] by Laplace transform (i.e., semi- 
group) techniques generalizing the Fourier transform approach of [Ze.l-2,LZe]; for 
discussion/geometric interpretation, see [Z.2,MZ.l-2]. In our stability analysis, we 
will use only a small part of the detailed information given in the proposition, 
namely L p — > L q estimates on the time-decaying portion H + S + R of the Green's 
function Q (see Lemma 4.1, below). However, the stationary portion E of the 
Green's function must be estimated accurately for an efficient stability analysis. 
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Lemma 3.3. Under the assumptions of Proposition 3.2, there holds, additionally, 
(3.27) U 2 (A°(x))- 1 H(x,t;y) = 0, 

where n 2 := ^jj j ^ denotes projection onto the final r coordinates of G. 

Proof. Equivalently, we must show that U 2 (A°)~ 1 7l* = for all 1 < j < J. 
This is most easily verified by the intrinsic property of C* and 1Z* (readily seen 
from our formulae) that they lie, respectively, in the left and right kernel of Bq. 
For, this gives 



(3.28) 



S 2fc*J:=(jJ °)(A°)-^- 

= bu 2 (A )- 1 n*, 

yielding the result by invertibility of 6, (see condition (A3)). 



Lemma 3.3 quantifies the observation that, in U = (u, v) f coordinates, the "par- 
abolic" variable v experiences smoothing under the evolution of (2.1), whereas the 
"hyperbolic" coordinate v does not. (Recall that U = (A°)~ 1 G). 

Section 4. Linearized stability. 

We now show, for the Lax case under consideration, that linearized orbital stabil- 
ity follows immediately from the pointwise bounds of Proposition 3.2, thus partially 
recovering the result of [MZ.2] that was stated in Proposition 3.1. This analysis mo- 
tivates the nonlinear argument to follow in Section 5. 

We carry out our analysis with respect to the conservative variable G = A (x)U, 
i.e. with respect to linearized equations (3.2). Similarly as in [Z.2,MZ.l-2], define 
the linear instantaneous projection: 



(4.1) 



/+oo 
E(x,t;y)G (y) dy 
-oo 

=: -5{t)G'{x), 



where G denotes the initial data for (2.1), and G = G{U{x)) as usual. The 
amplitude 5 may be expressed, alternatively, as 



/+oo 
e(y,t)G (y)dy, 
-oo 



where 



(4.2) 



E(x,t;y) =: G'(x)e(y,t), 
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i.e., 

(4.3) e(y,t) := £ f errfn f V -±M ) - errfn f V -^M ) ) l~ 

for y < 0, and symmetrically for y > 0. 
Then, the solution G of (2.1) satisfies 

/+oo 
(H + g)(x,t;0)G o (y) dy, 
-oo 

where 

(4.5) g-.= s + R 

is the regular part and H the singular part of the time-decaying portion of the 
Green's function Q. 

Lemma 4.1. For weak shock profiles of (1.1), under assumptions (Al)-( A3), (H0)- 
(H3), there hold: 

/+oo 
g(;t;y)f(y)dy\ LP <C(l + t)-^ 1 - 1 ^\f\ Lg , 
-oo 

/+°° _ , 
Q y (;t;y)f(y)dy\ LP < C(l + t)-^- 1 '^ 2 ^ + Ce~*\f\ LP , 
-oo 

/+oo 
GA-,t;y)f(y)dy\ LP <c(i + t)-^ 1 - 1 ^\f\ Lq + Ce-v t \f\ LP , 
-oo 

/+oo 
g xy (-, t- y)f(y)dy\ LP < C(l + f)-** 1 - 1 /')"^/^ + Ce^|/| r ,„ 
-oo 

and 

/+oo 
-oo 

for all t > 0, some C, r\ > 0, for any 1 < q < p (equivalently, 1 < r < p) and 
f E L q D W 1 *, where 1/r + l/q = l + 1/p. 

Proof. Bounds (4.6)-(4.9) follow by the Hausdorff- Young inequality together 
with bounds (3.16) and (3.17)-(3.20), precisely as in [Z.2,MZ.l-2]. Bound (4.10) 
follows by direct computation and the fact that particle paths Zj(y, t) satisfy uniform 
bounds 

l/C<\(d/dy) Zj \<C, 

for all y, t, by the fact that characteristic speeds aj(x) converge exponentially as 
x — > ±oo to constant states. B 
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Corollary 4.2. Let U be a weak shock profile of (1.1), under assumptions (Al)- 
(AS), (H0)-(H3). Then, strong spectral stability implies L 1 n LP — > LP linearized 
orbital stability, for any p > 1. More precisely, for initial data Uq G L 1 n L p , i/ie 
solution U = (u,v) t (x,t) of (2.1) satisfies the linear decay bounds 

(4.11) |Z7(-,t) + 5(t)^(-)U- < C(l + *)-* (1 " 1/p) (l^oki + |C/oj^)- 
Moreover, provided Uq G W 1,p , there hold also the derivative bounds 

(4.12) |(„(.,f) + 5 (fy(.))J LP <c(i + t)-^ 1 - 1 / p )(|t/oUi + l^oki.p). 

Proof. For (4.11), it is equivalent to show that, for initial data Gq G L 1 n L p , 
the solution G(x,i) of (3.2) satisfies 

(4.13) < C(l + t)-^ (1 - 1/p) (|G | L i + |Go| w i,„). 

But, this follows immediately from (4.4) and bounds (4.6) and (4.10), with q = p. 
Likewise, for (4.12), it is equivalent to show that, for initial data G G L 1 n W 1,p , 
the solution G(x,t) of (3.2) satisfies 

(4.14) Kna^ )- 1 ^-,*)-^-,*)))^^ <C(l + t)-^ 1 - 1 / p >(|GoUi + |Go| W i.p). 
This follows, similarly, from the derivative analog 

(4-15) 

(n 2 (A°)- 1 (G(x,t)-^(x,t))) a! = y °°(n 2 (A )- 1 (^ + ^)) ;i; (x,t;0)Go(2/)rf2/ 
of (4.4), together with (3.27) and bounds (4.6), (4.8), with q = p. B 

Section 5. Nonlinear stability. 

Finally, we establish our main result of nonlinear orbital stability with respect to 
perturbations Uq G L 1 nH 3 of weak viscous profiles (necessarily Lax type, by Propo- 
sition 1.1) of dissipative, symmetric hyperbolic-parabolic systems of type (1.1), of 
strength e > sufficiently small with respect to the parameters of the system in 
question. We follow the basic iteration scheme of [MZ.1,Z.2]; for precursors of this 
scheme, see [Go.2,K.l-2,LZ.l-2,ZH,HZ.l-2]. 

Define the nonlinear perturbation 

(5.1) U(x,t) :=U(x + 5(t),t) -U(x), 

where 5(t) (estimating shock location) is to be determined later; for definiteness, 
fix 5(0) = 0. Substituting (5.1) into (1.1), we obtain 

(5.2) G(U) t + F{U) X - (B(U)U X ) X = SG(U) X , 
and thereby the basic perturbation equation 

(5.3) (G(U)-G(U)) +(f(U)-F(U)) -(b(U)U x -B(U)U x ) = S(t)G(U) x , 

where U now denotes U(x + 5(t) : t) and U denotes U(x). 

5.1. Energy estimates. We begin by establishing the following basic energy 
estimate, by which we will eventually close our nonlinear iteration argument: 
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Proposition 5.1. Let Uq G H 3 , and suppose that, for < t < T , both the supre- 
mum of \8\ and the H 2 n W 2,oc norm of the solution U = (u,vY of (5.1)-(5.3) 
remain bounded by a sufficiently small constant £ > 0, for a small- amplitude shock 
profile as described in Proposition 1.1 of a system (1.1) satisfying (Al)-(A3), (H0)- 
(H3), with e := \U+ — U-\ sufficiently small. Then, there hold the bounds: 
(5-4) 

\U\ H *(t), J {\U X \ H 2 + \v x \ 2 H3 )(s)ds < C J (|C/| L oo(|t/| L oo + \U\ 2 L2 ) + 5 2 ) (s)ds, 
for allO<t<T. 

Remark 5.2. Note, with the expected decay rates 1^(^)1 ~ C(l + t) -1 / 2 , 

and \U\l2 ~ C(l +t -1 / 4 ), that the righthand side of (5.4) becomes order log(l +t), 
very nearly recovering the order one bound available in the constant-coefficient case 
[Kaw]. 

Lemma 5.3 ([SK]). Assuming (Al), condition (A2) is equivalent to either of: 
(Kl) There exists a smooth skew- symmetric matrix K{u) such that 

(5.5) Re (K(A°)- 1 A + B)(u)>e>0. 

A , A, B as in (2.2). 

(K2) For some 9 > 0, there holds 

(5.6) Re a(-^(A°)" 1 A( W ) - ^(A )- 1 B(u)) < -^| 2 /(1 + l£| 2 ), 
for all £ e R. 

Proof. These and other useful equivalent formulations are established in [SK] . B 

Lemma 5.4. Let (Al)-(A3) and (H1)-(H3) hold, and let U(x — st) be a viscous 
shock solution such that the profile {U(z)} lies entirely within a sufficently small 
neighborhood V C U of U*, and the speed s lies within a sufficiently small neigh- 
borhood of a p {U*): i.e., a profile as described in Proposition 1.1. Then, letting 
e := \U + — U-\ denote shock strength, we have for q = 0, . . .4 the uniform bounds: 

\dm(x)\ < Ce q+1 e- 0£ \ x \ 
(5 7) _ _ 

\d q x U\ L v < Ce q+1 ~ 1/p , 

for some C, 9 > 0. 

Proof. Though the bounds (5.7) are not explicitly stated in [Fre.2], they follow 
immediately from the detailed description of center manifold dynamics obtained in 
the proof, exactly as in the strictly parabolic case [MP] . B 
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Proof of Proposition 5.1. Writing (5.1) in quasilinear form 

(5.8) (A°U t - A°U t ) + (AU X - AU^ - [BU X - BU X ) = S(t)A°U X7 

where 
(5.9) 

1° := A°(U),A° := A°(U); A := A(U),A := A(U); B := B(U),B := B(U), 
using the quadratic Leibnitz relation 

(5.10) A 2 U 2 - A 1 U 1 = A 2 (U 2 - Ux) + (A 2 - A 1 )U U 

and recalling the block structure assumption (A3), we obtain the alternative per- 
turbation equation: 

(5.11) A°U t + AU X - (BU X ) X = M X {U)U X + (M 2 (U)U X ) X + S(t)A°U x + S(t)A°U x , 
where 

(5.12) M 1 (U):=A-A=^ dA(U + 6U)d9 S jU, 
and 

M 2 (U) :=B-B 
(5-13) f o \ 

\o (/o db(u + eu) de)u ) ' 

with 

(5.14) b:=b(U),b:=b(U). 

We now carry out a series of successively higher order energy estimates of the 
the type formalized by Kawashima [Ka.l]. The origin of this approach goes back 
to [K,MN] in the context of gas dynamics; see, e.g., [HoZ] for further discus- 
sion / references. 

Let K denote the skew-symmetric matrix described in Lemma 5.3 associated 
with A , A, B. Then, regarding A , K, we have the bounds 

, A° x = dA°(U)U x , K x = dK(U)U x , A x = dA(U)U x , B x = dB(U)U x , 

(5.15) _ 

A° = dA°(U)U t , K t = dK{U)U u A t = dA{U)U u B x = dB(U)U u 
and (from defining equations (5.1)-(5.2)): 

(5.16) \U X \ = \U X + U X \<\U X \ + \U X \ 
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and 

\U t \<C(\U x \ + \v xx \ + \5\\U x \) 
(5-17) < C(\U X \ + \U X \ + \v xx \ + \v xx \ + \S\\U X \ + \S\\U X \) 

< C(\U X \ + \U X \ + \v xx \ + \v xx \). 

Thus, in particular 

(5.18) \S\, \A% \K X \, \A X \, \B X \, \A% \K t \, \A t \, \B t \ < C(( + e). 

H 1 estimate. We first perform a standard, "Friedrichs-type" estimate for sym- 
metrizable hyperbolic systems. Taking the 1? inner product of U against (5.11), we 
obtain after rearrangement /integration by parts, and several applications of Young's 
inequality, the energy estimate 
(5.19) 

1 -(U,A U} t = (U,A U t }+ 1 -(U,A° t U} 

= -(U,AU X ) + (U, (BU X ) X ) + (U^M^U)^) + (U, (M 2 (U)U X ) X ) 
+ 5(t) (U, A°U X ) + 8(t) (U, A°U X ) + \ (U, A° t U) 

= \{U,A X U) - (U X ,BU X ) + (U,M 1 (U)U X ) - (U X ,M 2 (U)U X ) 

- U{t) (U, A° X U) + S(t) (U, A°U X ) + \ (U, A° t U) 
<-(U x ,BU x ) 

+ C J ((|^| + |^| + |^|)|t/| 2 + |^||t/||^| + |5||t/||^|) 

<-(U x ,BU x ) 

+ c J ((\u x \ 2 + \u\ 2 + \ Vx \ 2 + \v xx \ 2 )(\u\ + \u x \) + H 2 |tg) 

< -(U X ,BU X ) 

+ C (\U\ L ~{\U\ L ~ + \U\h) + |5| 2 ) + C(e + (\U X \ 2 L2 + \v xx \ 2 L2 ) , 

where e, ( > is as in the statement of the Proposition. Here, we have freely used 
the weak shock assumption and consequent bounds (5.7), as well as (5.15). (Note: 
we have also used in a crucial way the block-diagonal form of M 2 in estimating 
\(U X ,M 2 (U)U X )\ <Cf KH^H^xl in the first inequality). 

Likewise, differentiating (5.11), taking the L 2 inner product of U x against the 
resulting equation, and substituting the result into 

(5.20) l(U x ,A°U x ) t = (U x , (A°U t ) x ) - (U x ,A° x U t ) + ^(U x ,A° t U x ), 
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we obtain after an integration by parts: 
(5.21) 

^(U x ,A°U x ) t = -(U x , (AU X )) + (U, (BU X ) X ) + (U.M^U)^) + (U, (M 2 (U)U X ) X ) 

+ 5(t)(U,A U x } + 5(t)(U,A°U x } + ^{U,A° t U} 

= -\{U Xl A x U x ) - (U XX ,BU XX ) - (U XX ,B X U X ) + (U x , {M\{U)U X ) X ) 

- (U XX ,(M 2 (U)U X ) X ) + U(t)(U x ,A° x U x ) + d(t)(U x ,A U xx ) 

+ 5(t)(U x ,A° x U x ) - (U x ,A° x U t ) + ^(U x ,A° t U x ), 

which by (5.18), plus various applications of Young's inequality yields the next-order 

energy estimate: 

(5.22) 

^(U x ,A°U x ) t <-(U xx ,BU xx ) 

+ C J ((|t/| 2 + \d\ 2 )(\U xx \ + \U X \) + (e + ()(\U X \ 2 + \U x \\v xx \j) , 

< —(U xx , BU XX ) 

+ C (dC/lloo + |5| 2 ) + C(e + C) {\U X \ 2 L2 + \v xx \ 2 L2 ) 

1 ~ 

< — ^ \U X xi BU XX ) 



+ c([\U\ 2 Lac + \5\ 2 )+C{e + Q\U, 



2 

x\L 2 ' 



Next, we perform a nonstandard, "Kawashima-type" derivative estimate. Tak- 
ing the L 2 inner product of U x against K(A°)~ 1 times (5.11), and noting that 
(integrating by parts, and using skew-symmetry of K) 

^(u x , KU) t = ±(u x , ku t ) + \{u xtl ku) + ±(u x , k t u) 
= hu x ,ku t )-hu t ,ku x ) 

(5.23) 2 2 

--(U t ,K x U x ) + -(U x ,K t U) 

= (u x ,ku t ) + ^(u,k x u t ) + ^(u x ,k t u), 



we obtain by calculations similar to the above the auxiliary energy estimate: 
1 

(5.24) 2 



1 (U x ,kU) t <-(U x ,k(A°)- 1 AU, 



XI 

+ C(C + e)\U x \ 2 L2 + CC l \v xx \h + C(\U\l°o + \5(t)\ 2 ). 
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Adding (5.19), (5.24), and (5.22) times a suitably large constant M > 0, and 

recalling (5.5), we obtain, finally: 

(5.25) 

X - (V, A°U)+(U X , KIT) + M(U X , A°U x )) t 

< -0(\U x \l* + \v xx \\ 2 ) + C (\U\ L ~{\U\ L ~ + \U\ 2 L2 ) + |5| 2 ) , 

9 > 0, for any £, e sufficiently small, and M, C > sufficiently large. 

Higher order estimates. Performing the same procedure on the once- and twice- 
differentiated versions of equation (5.11), we obtain, likewise, the H q estimates, 
q = 2, 3, of: 

^(dr'uA^r'u) + (d^ 1 u x ,kdr 1 u) + M(d«u,A°d«u)) t 

(5.26) <-0(mU\l 2 + 18^11.) 

+ (e + C)\U x \ 2 Hq -2 + C (|tf| L oo(|tf| L oo + \U\ 2 L2 ) + |5| 2 ) . 

We omit the calculations, which are entirely similar to those carried out already. 

Final estimate. Summing our H q estimates from q = 1 to 3, and telescoping the 
sum of the righthand sides, we thus obtain, for e, £ sufficiently small: 

3 

(5.27) g =i Z 

< -9(\U X \ H 2 + \v x \ 2 H3 ) + C + \U\ 2 L2 ) + |5| 2 ) , 

or, integrating from to t: 

3 



J2 {(A°dr 1 u,dr 1 u)+hd*u,kd*- 1 u) + M{A d^d q u))\i 

(5 - 28) <- f 9{\U x \ H , + \v x \ 2 m ){s)ds 

Jo 

+ C J (\U\ LO o(\U\ L <*> + \U\ 2 L2 ) + |5| 2 ) (s)ds. 

Noting that, for M sufficiently large, we have for each q, by Young's inequality, 
and positive definiteness of A : 
(5.29) 
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for some 6 > 0, we may rearrange (5.28) to obtain our ultimate goal: 

\U\ 2 H3 (t)+ I 0{\U x \ H 2 + \v x \ 2 H3 ){s)ds 

(5.30) Jo t 

< C\U\ 2 H3 (0) + C (\U\ L °o{\U\ L ~ + \U\ 2 L 2) + |5| 2 ) (s)ds, 

from which the result immediately follows. | 

5.2. Nonlinear iteration. We now carry out the nonlinear iteration, following 
[Z.2,MZ.l]. For this stage of the argument, it will be convenient to work again with 
the conservative variable 

(5.31) G := G(U) - G(U), 
writing (5.1) in the more standard form: 

(5.32) G t -LG = N(G, G x ) x + 6(t)(G x + G x ), 

G := G(U), where 

N(G,G x ) = 0(\G\ 2 + \G\\v x \), 
N(G,G X ) X = 0(\G\ 2 + \G\\v x \ + \G x \\v x \ + \G\\v xx \), 

so long as |67|, \G X \ remain bounded. Here, v denotes the second coordinate of the 
alternative perturbation variable U = (u,vY defined in (5.1). 
By Duhamel's principle, and the fact that 

/oo 
G(x,t; y)G x (y)dy = e Lt G x (x) = G x (x), 
-00 

we have, formally, 

/oo 
g(x,t;y)G (y)dy 
-00 



(5.35) 



pt poo 

+ / / g y (x,t-s;y)(N(G,G x ) + 6G)(y,s)dyds 

JO J -00 



'0 J -00 

+ 6(t)G x . 



Defining, by analogy with the linear case, the nonlinear instantaneous projection: 
<p(x,t) :=S(t)G x 

) 

E(x,t;y)G (y)dy 



(5.36) 



00 

/ E y (x, t - s; y)(N(G, G x ) + SG)(y, s)dy, 
^-00 
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or equivalently, the instantaneous shock location: 

/oo 
e(y,t)Go(y)dy 
-oo 



(5.37) - x 



ft f-\-oo 

+ / / ey(y,t-s)(N(G,G x ) + 6G)(y,s)dyds, 

JO J-oo 



where E, e are defined as in (3.15), (4.3), and recalling (4.5), we thus obtain the 
reduced equations: 

/oo 
(H + g)(x,t;y)G (y)dy 
-oo 

ft fOO 

(5.38) +/ / H(x,t-s;y)(N(G,G x ) x + 8G x )(y,s)dyds 

JO J-oo 



t poo 

G y (x, t-s-y) (N(G, G x ) + SG) (y, s)dy ds, 



'0 J-oo 

and, differentiating (5.37) with respect to t, 



(5.39) 



/oo 
e t (y,t)G (y)dy 
-oo 
pt f-\-oo 

+ e yt (y,t- s)(N(G,G x )+5G)(y,s)dyds. 

Jo J-oo 



Note: In deriving (5.39), we have used the fact that e y (y, s) —r as s — > 0, as 
the difference of approaching heat kernels, in evaluating the boundary term 



/-t-oo 
e y (y,0)(N(G,G x ) + 5G)(y,t)dy = 0. 
-oo 

(Indeed, \e y (-, s)\ L i — > 0, see Remark 2.6, below). 

The defining relation S(t)u x := —ip in (5.36) can be motivated heuristically by 

g(x, t) - <p(x, t)~G=(^y x + s( t ) ,t)-(*y X ) 

~G{x,t) + 6{t)G x {x), 

where G denotes the solution of the linearized perturbation equations, and G the 
background profile. Alternatively, it can be thought of as the requirement that the 
instantaneous projection of the shifted (nonlinear) perturbation variable G be zero, 
[HZ. 1-2]. 
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Lemma 5.5 [Z.2]. The kernel e satisfies 

(5.41) \e y {-,t)\ LP ,\e t {;t)\ LP <Cr^ 1 - 1 ^\ 

(5.42) |e tv (-,t)| L p<Cr'( 1 - 1 / p )- 1 / 2 , 

/or all t > 0. Moreover, for y < we have the pointwise bounds 

(5.43) \e y (y,t)\,\e t (y,t)\<Ct-h- 



! (y+g_t) 2 



Mi 



(y+a_t) 2 



(5.44) |e^(y,t)| < Ct"^ 

for M > sufficiently large (i.e. > 4b±), and symmetrically for y > 0. 
Proof. For definiteness, take y < 0. Then, (4.3) gives 

(5.45) e tf (y, £) = ( * J (K(y + a.t, t) - K(y - a.t, t)) , 

\ u + U— / 

(5.46) e t (y, t) = ( i ) ((K + K y )(y + a_t, t) - (K + K y )(y - a_t, t)) , 

(5.47) 

e* y (y,*) = ( ) ((K y + K yy )(y + a-t,t) - (K y + K yy )(y - a-t,t)) , 

where 

-y 2 /4b_t 

(5.48) K(y,t) 



denotes an appropriate heat kernel. The pointwise bounds (5.43)-(5.44) follow 
immediately for t > 1 by properties of the heat kernel, in turn yielding (5.41)- 
(5.42) in this case. The bounds for small time t < 1 follow from estimates 



\Ky(y + a-t,t)-K y (y-a-t,t)\ = \ / K yy (z,t)dz\ 

J y+a-t 

( 5 - 49 ) <Ct~^ 2 e^dz 

J v+a-t 



'y+a-t 

< Ct~ 1 / 2 e- : 
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and, similarly, 

/a_ t 
Kyyy(Z,t)dz\ 
-a-t 



y-a-t 



(5.50) < ct - 2 r ~\^ dZj 

Jy+a-t 



(y+a_t) 2 



< Ct~ x e mj 

The bounds for \e y \ are again immediate. Note that we have taken crucial account 
of cancellation in the small time estimates of et, ety D 

Remark 5.6: For t < 1, a calculation analogous to that of (5.49) yields \e y (y, t)\ < 

(y+a,_t) 2 

Ce Mt , and thus |e(-, s)\ L i — > as s — > 0. 

With these preparations, we are ready to carry out our analysis: 
Proof of Theorem 1.3. Define 

t(t) := SUp [ (\U(; S)\ LP + \V X (; s)\ LP )(l + fl)** 1 " 1 /*) 

(5 51) 0<s<t,2<p<oo L 

+ \5(S)\(1 + S) 1 ' 2 + {\5(S)\ + \U(; S)\ H , + \U(; S )\ W ^) ] . 

We shall establish: 

Claim. For all t > for which a solution exists with £ uniformly bounded by 
some fixed, sufficiently small constant, there holds 

(5-52) CW<C 2 (Co + CW 2 ). 

From this result, it follows by continuous induction that, provided Co < 1/4C 2 , 
there holds 

(5.53) C(t) < 2C 2 Co 

for all £ > such that £ remains small. By standard short-time theory /continuation, 
we find that the solution (unique, in this regularity class) in fact remains in H 2 for 
all t > 0, with bound (5.53), at once yielding existence and the claimed sharp L p 
bounds, 2 < p < oo. Thus, it remains only to establish the claim above. 



\v x 
by 



Proof of Claim. We must show that each of the quantities |{7|lp(1 + s) 2 ^ 1 l / p \ 
c Ml + s)s (1 - 1/p) , I^Kl + s) 1 / 2 , \5\, \U\ H 2, and \U\ ^2,00 are separately bounded 



(5.54) 



^(Co + CW 2 ), 
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for some C > 0, all < s < t, so long as ( remains sufficiently small. 

Provided that we can establish the others, the final two bounds follow easily from 
the energy estimates of Proposition 5.1. For, by bounds (5.54) on \S\ and \U\loo, 
and the assumption that |L/|#-2 nW /2,oo < ( remains small, we obtain from (5.4) that 

(5.55) \U\m(s) < C 3 (Co + C(t) 2 ) log(l + S), 

for < s < t. Interpolating with assumed bound 

l^|L 2 («)<^ 3 (Co + C(t) 2 )(l + «)- 1/4 , 



we obtain 



(5.56) 



2^irr|2/6| rr |l/6 



\u\ H ,(s)<c 3 (Co + atr)\u\^\u 



L 



2 



<c 3 (Co + CW 2 )(i + t)- 1/25 , 

and, by Sobolev estimate, 

(5.57) \U\ W 2^ < \U\)g\U\)g < C 3 (Co + C(t) 2 )(l + *)" 1/5 \ 

both better than claimed. 

Thus, in order to establish the result, we have only to establish the remaining 
bounds, on |C/|x, P , |i> x |lp, and \S\. These will be carried out using the Green's 
function estimates of Lemma 4.1. Accordingly, we first convert the problem to 
conservative, G coordinates, via: 

Observation 5.7. It is sufficient to establish corresponding bounds on 
\ v g,x\lp, 1 5 1, and 1 5 1, where vq(x, t) := A°(x)~ 1 G(x, t). 

Proof. We have 

U = G-\G) - G~\G) = (A°)zX(x, t)G 

(5.58) / r 1 \ 

:= U (A°)- 1 (U(x) + 6U(x,t))de) G, 

where 

(5.59) l(^ )ave(M) - A ^)" 1 ! <C\U\< C\G\ 
and 

(5-60) | {(A°):; e (x,t) - A\x)~ 1 ) x I < C(\U X \ + \U X \) < (, 



whence 



U lp < C \G\i,v 
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and 

\U X - ((A°)~ 1 G) X \ LP < C(\G\ L ,\G X \ L - + C\G\ LP ) < CC\G\ LP , 
from which the result easily follows. B 

By (5.38)-(5.39), we have 

/oo 
(H + g)(x,t;y)G (y)dy\ LP 
-oo 
rt p oo 

+ 1/ / if(a:,t- S ;y)(iV(G,G x ) x + (jG' ;E )( 2 /, S )rf 2 /^| L P 
(5.61) Jo i-oo 

/•t />oo 

+ 1 / / ^ y (a;,t-s;y)(iV(G',G a; ) + 5G)(y,s)^rfs| LP 

JO J-oo 
=: I a + h + Ic, 

(5.62) 

\v Gl x\L*(t) :=\{(A )- 1 G)J LP (t)<\ J°° ( K (A°)- 1 g')jx,t;y)Go(y)dy\ LP 



•>t poo 

+ ' 



/ / ((A°)- 1 g y ) x (x,t-s;y)(N(G,G x ) + SG)(y,s)dyds\ LP 

JO J-oo 



=: Ha + Ih, 



(5.63) 



oo 

1*1(0 < I / e t (y,t)G (y)dy\ 

-oo 
ft p+oo 



pi p-too 

+ | / / e yt (y,t-s)(N(G,G x ) + 6G)(y,s)dyds\ 

Jo J-oo 



and 



(5.64) 



'0 J-oo 
=: III a + III b , 



\S\(t)<\ I e(y,t)G (y)dy\ 

-oo 

»t p+oo 

+ ' 



pt r+oo 

/ / e y (y,t-s)(N(G,G x ) + 8G)(y,s)dyds\ 

Jo J-oo 



=: IV a + IV b . 

We estimate each term in turn, following the approach of [Z.2,MZ.l]. The linear 
term I a satisfies bound 

(5.65) I a <CCo{l + t)-^ l - l ^\ 
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as already shown in the proof of Corollary 4.2. Likewise, applying the bounds of 
Lemma 4.1 together with (5.33) and definition (5.51), we have 



(5.66) 



h = \ / H(x,t-s;y)(N(G,G x ) x + 8G x )(y,s)dyds\ LP 

JO J-oc 

< C I e-^- a \\G\ Loo + \v G , x \ L °° + \S\)\G\ W 2, P (s)ds 

Jo 

< C((t) 2 [ e-^ t - s \l + s)- 1 ' 2 ds 

Jo 

<C7C(t) 2 (l+t)" 1/2 , 



and (taking q = 2 in (4.7)) 



ft />oo 

/ G y (x,t- s;y)(N(G,G x ) +5G)(y,s)dyds\ LP 

J — oo 

< C f e-^- s \\G\ L ~ + \v G , x \ L ~ + \8\)\G\ LP (s)ds 
Jo 

+ C f (t- s)- 3 / 4+1 / 2p (\G\ L ~ + \v G , x \ L ~ + \5\)\G\ L2 (s)ds 
Jo 

<C((t) 2 f e -v(ts)( 1 + s )-h(i-i/p)-i/2 ds 
Jo 

+ C((t) 2 [\t - S )- 3 / 4+1 /^(l + s)- 3 /*ds 
Jo 



(5.67) 



|(i-Vp) 



<caty(i+t)-z 



Summing bounds (5.65)-(5.67), we obtain (5.54), as claimed, for 2 < p < oo. The 
desired bounds on II a and lib follow by an identical calculation, once we notice 

that, in Lemma 4.1, ((A )' 1 ^ and ((A )' 1 ^ satisfy the same L q -> U> bounds 

as do (H + Q) and (H + Q) y , respectively. 

Similarly, applying the bounds of Lemma 5.5 together with definition (5.51), we 
find that 

/oo 
e t (y,t)G (y)dy\ 
-oo 

(5 - 68) <|et(y,f)U~(f)|GoUi 

<CCo(l + t)" 1/2 
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and 



(5.69) 



while 



and 



(5.71) 



IIh = \ / e yt (y,t-s)(N(G,G x ) + 5G)(y,s)dyds\ 

JO J-oc 

< / \e y t\Li(t - s)(\G\ L oo + \v G , x \ L °o + \S\)\G\ L 2(s)ds 
Jo 

<C((t) 2 [\t- s )- 3 /\l + s)- s / 4 ds 
Jo 

<cc(t) 2 (i + t)" 1/2 , 



/oo 
e(y,t)G (y)dy\ 
. ... _ -oo 

{ ] < Hy,t)\ L At)\GoW 

<CCo 



IV b =\ / e y (y,t - s)SG(y, s)dyds\ 

Jo J -oo 

< I \e y \ L ^(t - s)(\G\ L °° + \vg,x\l°° + \S\)\G\ L 2(s)ds 
Jo 



< cat) 2 [ (t-s)- i/4 (i + s )- s / 4 ds 

Jo 

< cat) 2 . 

Summing (5.68)-(5.69) and (5.70)-(5.71), we obtain (5.54) as claimed. 
This completes the proof of the claim, and the result. | 
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